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TECHNICAL MEMORANDUM X-64534 


INTERPOLATION OF BODE'S TABLES OF THE 
SOLAR CONTINUOUS ABSORPTION COEFFICIENT 

SUMMARY 


Knowledge of the total continuous absorption coefficient for a wide 
range of wavelengths is vital to many areas of solar physics, such as the 
construction of model atmospheres and line profile calculations. Recently, 
Gerhard Bode [1] compiled a set of tables for the continuous absorption 
coefficient for the sun covering a wide range of wavelengths , electron 
pressures, and temperatures. To obtain the absorption coefficient as a 
function of electron pressure and temperature at an untabulated wavelength, 
one can employ standard interpolation formulas using the several wavelengths 
listed, which are between major absorption edges and span the desired wave- 
length. At a given wavelength the absorption coefficients at particular 
temperatures and pressures not given in the table can be obtained by a double 
linear interpolation scheme. This paper describes a computer subprogram 
which has been developed to calculate the solar continuous absorption 
coefficient for any given model atmosphere. 


INTRODUCTION 


The subprogram used to interpolate the Bode data employs a method 
of two-way differences. This method is described in detail by Scarborough 
[ 2 ]. 


It is hoped that an example problem of the technique of two-way dif- 
ferences will clarify the procedure used to interpolate the Bode tables by 
computer. A flow chart (Appendix A) is also used to indicate the processes 
undertaken by the computer; the subroutine is listed in Appendix B. 



THE TECHNIQUE OF DOUBLE INTERPOLATION 


The example problem uses a table of values of log (electron pressure/ 
absorption coefficient) at 5250. 20 A vs log (electron pressure) and 0 , 
where 9 is 5040/absolute temperature. 

It should be noted that there is no Bode table for 5250. 20 A; rather 
Table 1 was com|)iled using a simple straight-line interpolation of a Bode 
table at 4500. 00 A and one at 5698. 50 A 


TABLE 1. VALUES OF LOG (Pe/K) 



Log (Pe) 

0 

-2. 5 

-2. 0 

-l.,5 

-1. 0 

-0,5 

0. 0 

1. 30 

24. 376 

24. 569 

24. 658 

24. 693 

2A713 

24. 758 

1. 25 

24. 412 

24. 625 

24. 727 

24. 767 

24. 784 

24. 806 

i. 20 

24. 443 

24. 681 

24. 797 

24. 884 

24. 861 

24. 872 

1. 15 

24. 465 

24. 729 

24. 865 

24. 922 

24. 942 

24. 950 

1. 10 

24. 452 

24.767 

24. 927 

24. 999 

25. 024 

24. 033 

1. 05 

24. 313 

24. 768 

24. 985 

25. 072 

25. 105 

25. 117 

1. 00 

23. 914 

24. 627 

25. 004 

25. 138 

25. 183 

25. 202 

0. 95 

23. 298 

24. 200 

24. 879 

25. 174 

25. 255 

25. 281 


Suppose that a value of log { Pe/K) is needed for particular values of 
9 and log (Pe) not given in Table 1. The method of two-way differences 
would then be used as follows: 

Problem: Log (Pe/K) needs to be calculated at a 0 value of i. 22 
and a log ( Pe) value of -2. 2 at 5250. 20 A 

Solution: Tables 2A and 2B show forward differences for constant 
log (Pe) values, and Tables 3A and SB show differences for constant 0 
values (hence, two-way differences) . A general formula for double inter- 
polation uses these differences in providing the desired value of log (Pe/K) . 
Note that log (Pe/K) = 24. 412 at 0 = 1. 25 and log (Pe) = -2, 5 is used as a 
"starting point" for taking differences. The starting point is that value of 
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log (Pe/K) in Table i for which the corresponding values of log (Pe) and 9 
are closest to and above those which are being interpolated. 


TABLE 2A. CONSTANT LOG (Pe) VALUES; LOG (Pe) = -2. 5 



K(0, -2.5) 

1+0® 

* “eo 

2+0^ 

* “eo 

6 Q = 1* 25 

24. 412 

0. 031 


6i = 1. 20 

24. 443 

0. 022 

-0. 009 

02 = 1. 15 

24. 465 




TABLE 2B, CONSTANT LOG (Pe) VALUES; LOG (ge) = -2. 0 



K(0, -2.0) 

l+o“ 

* “ei 

. 2 + 0 “ 

00 = 1. 25 

24. 625 

0. 074 


01 = 1. 20 

24. 699 

0. 030 

-0. 044 

02 = 1. 15 

24. 729 




a. Kqo is the starting point for taking differences; 
K(0, Log Pe) - Log (Pe/K) 


A K 


ePe 


= First-order difference, 


.(h-i 
A K 


0Pe 


A 2+0 

A K 


A 0+2 

A K 


OPe 


0Pe } 


Second-order difference. 


Also, 


, 1+1 , 1+0 , 1+0 , 0+1 , 0+1 

A = A . - A K„= A K._ - A K 

0Pe 01 60 iPe OPe 
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TABLE 3A. CONSTANT 6 VALUES; 6 » L 25 



K( 1. 25, Log Pe) 

. 0 + 1 “ 

.0+2^ 

'^IPe 

Log Peg = -2. 5 

24. 412 

0. 213 


Log Pe^ = -2. 0 

24. 625 

0. 102 

-0. Ill 

Log Pe 2 = -1. 5 

24. 727 




TABLE 3B. 

CONSTANT 0 

VALUES; 0 

- 1. 20 


KU.20, Pe) 

. 0 + 1 “ 

^ “iPe 

0+2^ 

IPe 

Log Peg = -2. 5 

24. 443 

0. 256 


Log Pej = -2. 0 

24. 699 

0. 098 

-0. 158 

Log Pe 2 = -1. 5 

24. 797 




Kqo is the starting point for taking differences; 
K(e, Log Pe) = Log (Pe/K) 


A K 


0 Pe 


= First-order difference, 


A K 


0Pe 


. 2+0 
A K 


0Pe 


= Second-order difference. 


0+2 
A K 


0Pe 


Also, 

1+1 ^ 1+0 ^ 1+0 . 0+1 ^ 0+1 
A =A «■ -A K = A K,^ -A K 

0Pe 01 00 IPe OPe 


The general formula for double interpolation using two-way differences 




















• where 


K = Koo + u • 


A 

A Kqo + V 


A 


0+1 


^00 


+ 


2 ! 


u(u-l) 


+ 2u • V 



+ v(v-l) 



6 - e 


u 


1 . 


01 - 00 


and 


V = 


Log (Pe) - Log (Pe)p 
Log ( Pe; 1 - Log ( Pe) q 

From Table 1, 


Koo = 24.412 

A Koo = 

0. 031 

A Koo - 

0. 213 

A 2+0 

A Koo - 

-0. 009 

0+2 

A Koo - 

-0. Ill 


Also, 


i + i A A A AAO 

A Koo = A - A k:^„= 0.043 


'01 


00 


2+0 

• A Koo 



For u and v : 


1. 22 - 1. 25 „ 

“ = -0. 05 = 


V = , , , 

0 . 0 


Then, 


Log (Pe/K) 


= 24. 412 + ( 0. 6) ( 0. 031) + ( 0. 4) 
(0. 6) (-0. 4) (-0. 009) + 2( 
+ (0. 4) (-0. 06) (-0. Ill) 


1 

+ — 
2 


(0. 213) 

0. 6) (0. 4) (0. 043) 
= 24. 541 


Before the computer takes differences from Table 1 , it first finds a 
suitable starting point. In the preceding example the computer would go to 
24. 412 to start taking differences. If a 0 value of 1. 12 and a log (Pe) value 
of -1. 1 were read into the computer, the calculations would begin at 24. 865 
as a starting point. The computer scans those values of 6 and log { Pe) in 
Table 1 until it finds the ones that are closest to but higher in the table than 
the 6 and log ( Pe) values read in. Then the log ( Pe/K) value corresponding 
to these values of 6 and log (Pe) is used as a starting point for taking 
differences. In this way the best possible interpolated value of log (Pe/K) 
is obtained. 


CONCLUSIONS 


Some restrictions on using the subroutine should be mentioned. In 
Table 1 , the values of 9 must increase by equal amounts ( -0. 05 in this 
case) ; similarly, the values of log (Pe) must increase by equal steps. 

For the subroutine to function, it must be able to take first- and 
second-order differences from Table 1. As long as the value of 9 in the 
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interpolation is not equal to or smaller than 1. 00 and the value of log (Pe) 
is not equal to or greater than -0, 05, the calculations in the subroutine will 
be carried out. If this restriction is not observed, the table must be made 
larger, or the interpolation subroutine must be made to take only first- 
order differences. 
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CHART 


COMPILE VALUES OF 
THETA AND LOG (PE) 
FOR TABLE 


Actually, the values 
of log (Pe/K) are 
read into the main 
program; then, a 
COMMON statement 
is used. 











&FOR.IS OPAC 

SUBROUTINE OP AC ( X « THET A . PE . WV. AK.CG ) 

01 MENS T ON XKAPP r 9 .TrvXt,0GT7T VTHET 19) 

COMMON XKAPP 

THE! (1) =l.i)b 

90 DO ICO 1=2.9 

IOC THET( n=THET( I-l)-.05 

CONTINUE 

XLOGtr)=-T.O ■ " ^ — 

120 DO 130 K=2.7 

n-ST) XL0G(K)=XLUG(K-1 )+.b 

CONTINUE 

160 XTHET = 0.0 - ■ 

170 YTHET=0.0 

180 1 = 1 - - - -- 

190 1=1+1 

-?CrO XTHET - T H ET t I ) 

YTHET=THETA 

220 IF ( .05-(XTHET-YTHET)i -19O.?3O*23^0 

230 U=0.0 

240 U=(XTHET-YTMET)/.05 

250 K=1 

- 26 re 

270 ZL06=0.0 

280 ZLOG=XLOG(K) - 

290 WLOG=0.0 

WLOG= .434294+AL'0G(-P£) 

310 IF (.5+(ZL0G-WL0G) ) 260.320.320 

-5^0 v=^T-e 

330 V=(WL0G-ZL0G)/.5 

340 ■ FIRT=0.0 - — - - 

350 FIRT=XKAPP( I+1.K)-XKAPP( I.K) 

360 SECT=0.0 - 

370 SECT=XKAPP( 1+2 .K ) -XKAPP ( I+l.K) 

-3^^^ F ' l -R t -= 0 -re 

390 FIRL=XKAPP( I .K+1 )-XKAPP( I .K) 

400 SECL=0.0 

410 SECL= XKAPP ( I . K+2 ) -XKAPP ( I .K+1) 

420 SECTO = 0.0 - 

430 SECTO=SECT-FIRT 

-44^0 5 ECL 0 =0. 0 

450 SECLO=SECL-FIRL 

460 DA=0.0 - - - - 

470 DA = XKAPP ( I + 1 . K+1 ) -XKAPP ( I .K+1) 

480 DX=0.0 -- - - - 

490 DX=DA-FIRT 

-5313- X C AP=0.0 

510 XCAP=XKAPP( I .K)+U»FIRT+V*FIRL+(SECT0»(U-1. ) *U+2 .*U*V*DX+SECLO* ( V-1 
l.)*V)/ 2. 

XAKCG=0.0 

XAKCG=23.6203-(XCAP-WL06) - 

Y = 0.0 

Y i-xaK C G » 2. 30258 

AKCG=0.0 

10 AKCG=EXP(Y) 

RETURN 

END - ----- - ----- — 


COMPUTER SUBPROGRAM 
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